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Abstract: We study excitation spectra of BPS-saturated topological solutions — the 
kinks — of the scalar field model in (1 + 1) dimensions, for three different choices of the 
model parameters. We demonstrate that some of these kinks have a vibrational mode, apart 
from the trivial zero (translational) excitation. One of the considered kinks is shown to have 
three vibrational modes. We perform a numerical calculation of the kink-kink scattering in 
one of the considered variants of the model, and find the critical collision velocity Vcr 
that separates the different collision regimes: inelastic bounce of the kinks at Ujn > Ucr, and 
capture at v\n < Ter- We also observe escape windows at some values of fin < v^r where 
the kinks escape to infinity after bouncing off each other two or more times. We analyse 
the features of these windows and discuss their relation to the resonant energy exchange 
between the translational and the vibrational excitations of the colliding kinks. 
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1 Introduction and motivation 

Topological defects in (1+1) space-time dimensions have been a subject of active research [1- 
3]. Properties of these defects, their interaction with other defects and impurities are 
widely used to model various phenomena in condensed matter physics [4]. For instance, 
field models considering a real scalar field with polynomial self-interaction are routinely 
employed to model phase transitions. The most widely known example of such a model is 
probably the field theory. It has a topological solution — the kink — that can be used 
to describe a system undergoing a second order phase transition. More complex processes, 
when a series of consecutive phase transitions is to be modelled, make necessary the use of 
polynomial self-interactions with the degree higher than four. 

Field models that can develop topological solutions are also very important in the 
context of classical and quantum field theory, high-energy physics, cosmology, as well as 
hadron and nuclear physics [5-10]. In this relation, one has also to mention the progress in 
the research of properties of strings, vortices and monopoles [1, 11-23]. 

To stress the importance of studying topological defects in (1 + 1) dimensions, we 
note that, while being sufficiently easily treatable, they can correctly describe structures in 
higher dimensions, such as, for instance, a smooth domain wall in (3 + 1) dimensions whose 
profile can be described by a (1 + l)-dimensional kink. 

The existence of eigenmodes in the excitation spectrum of a solitary wave is related to 
its stability against small perturbations; these modes can also affect the interaction of the 
solitary waves with one another and with external objects. A soliton in a translationally 
invariant model can always be shifted by a constant vector, which means that there always 
exists the zero (translational) mode. Excitations with higher energies — vibrational modes 
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— can give rise to a rich collection of resonance phenomena in kink-(anti)kink collisions, as 
well as in the scattering of kinks off impurities [24-34] . 

This article deals with topological solitons of the (1 + l)-dimensional held the¬ 
ory. This model has been employed, in particular, to model massless mesons with self¬ 
interaction [35], and to describe isostructural phase transitions [36]. We study for the hrst 
time the excitation spectra of the kinks that occur in this model, for three different choices 
of the model’s self-interaction. We show that some of the studied kinks have a vibrational 
excitation mode. To illustrate how the excitation spectrum of a single kink affects the 
latter’s interaction with other kinks (and other spatial defects), we select one of the kinks 
that have a vibrational excitation and study its collisions with the corresponding antikink. 

Our study is structured as follows. Section 2 provides a brief introduction into general 
properties of static solutions with hnite energy in (1 -|- 1) dimensions. In Section 3, we 
deal with the kinks of the model (with three different choices of the model parameters) 
and perform a study of their excitation spectra. Section 4 presents the results of our 
numerical study of collisions between one of the kinks that are considered in Section 3 
and the corresponding antikink, and a discussion of the observed resonance phenomena. 
Section 5 delivers an outlook and the conclusion. 


2 Kinks in (1 + 1) dimensions 


We consider a held-theory system with a single real scalar held (p{t, x) in one spatial and 
one temporal dimensions, with the Lagrangian 

^ = \ /^ = 0 , 1 , ( 2 . 1 ) 

where the self-interaction of the held (/?, +(</?), is assumed to be bounded from below 
and hence can be thought of as a non-negative function of (p. The energy functional, 
corresponding to Eq. (2.1), is 



( 2 . 2 ) 


The Lagrangian (2.1) yields the following equation of motion for the held (p(t,x): 


^ dV ^ 
□ — 0 , 
dp> 


(2.3) 


where □ = is the d’Alembertian. For a static conhguration, (p = p>{x), this becomes 


(Pip _ dV 
dx"^ dip ’ 


(2.4) 


which can be transformed into 


dip 

dx 


^V^viip) 


(2.5) 





or further into 


( 2 . 6 ) 


_ ^dW 
dx dip ’ 

where the superpotential W(p) is related with V(p) via 



If the potential V(p) has two or more minima pi, P 2 , ■ ■ ■, yielding the same minimal 
value V(p^) = 0, i = 1,2,..., the energy of a static solution p = p{x) can be written as 


OO 

E = E'bps + 2 / 

— OO 


dp dW 
dx dp 


2 

dx. 


( 2 . 8 ) 


where 

^^BPS = \W[p{+ix)] - W[p{-oo)]\ (2.9) 

is the energy of the BPS-saturated solution [2, 37, 38], where the static field p(x) fulfils 
Eq. (2.6) and therefore the integrand in Eq. (2.8) turns into zero. The BPS-saturated solu¬ 
tion therefore has the smallest energy among all the static solutions interpolating between 
two given adjacent minima of the potential. 


p{—oo) = lim p{x), p{+oo) = lim p{x). 

x^—oo a:i^H-oo 

Note that the static solution p{x) has to converge sufficiently fast to one of the minima of 
the potential in order that the energy be finite. 


lim p{x) = Pi, lim p{x) = p^ (2.10) 

(this is also true if the field depends on time). As usual, solutions where p^ = Pj are 
called non-topological, whereas p^ ^ pj corresponds to a topological solution. We will call 
the family of all solutions with identical spatial asymptotics a “topological sector”, with 
the corresponding notation, e.g., {pi,pj) for the topological sector with the asymptotics of 
Eq. (2.10). 

The field theory being Lorentz-invariant, a static solution of Eq. (2.4) or Eq. (2.6) can 
be boosted to produce a soliton moving with a constant velocity v, with the energy given 

by 

\/l — 

with M being the energy of the static kink. We reserve the term “kink” (“antikink”) for 
topological BPS-saturated solutions that connect any two adjacent minima of the self¬ 
interaction potential, and their boosts. 

In order to analyse the excitation spectrum of a static kink p]^[x), we add to it a small 
perturbation, 

p{t,x) = p]^{x) + 5p{t,x), <C ||(/3k||, 
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( 2 . 11 ) 




( 2 . 12 ) 


and, taking in the equation of motion (2.3) terms linear in 6(p{t,x), obtain: 

■ 6(p = 0. 


d^6if d^Scp ^ (PV 




dP dx^ d(p^ 

Looking for x) in the form 

5ip{t^ x) = il^{x) cos{ujt), 

we obtain from (2.12) a boundary value problem 


d^V 
dx"^ d(p‘^ 




V’(x) = uP^{x), 


(2.13) 


where 'ip{x) has to satisfy the usual Schrodinger-like conditions, i.e., it has to be smoothly 
differentiable and square integrable over the real axis. The similarity to the Schrodinger 
equation can be further exploited by denoting 

d'^V 


U{x) = 


dip^ 




(2.14) 


making Eq. (2.13) an eigenvalue problem for the Hamiltonian 

H = -^ + U{x). 

It can easily be shown that the Hamiltonian (2.14) always has a zero eigenvalue, corre¬ 
sponding to the translational excitation mode. To demonstrate this, take the derivative of 
Eq. (2.4) with respect to x, which gives 


d^p' d?V I 


(2.15) 


If we substitute p = p]:^{x)^ this equation becomes an identity (as the kink is a solution of 
the equation of motion), at the same time coinciding with Eq. (2.13) if one selects a; = 0. 
Hence, 

?/)0(x) = <y9k(x) (2.16) 

is the eigenfunction of the Hamiltonian (2.14), corresponding to the eigenvalue w = 0. The 
square integrability of 'ipoix) follows from the fact that the energy of the kink is finite, cf. 
Eq. (2.2), therefore the eigenvalue w = 0 belongs to the discrete part of the excitation 
spectrum. 


3 The kinks of the model and their excitation spectra 

The p^ model is described by the Lagrangian (2.1), where the potential is a polynomial 
having the degree eight. The shape of the potential can vary depending on the values of 
the polynomial coefficients. In turn, different shapes of the potential can generate different 
sequences of phase transitions in condensed matter systems, as shown in ref. [39]. Our 
choice of the self-interaction parameters follows this reference; the specific potentials that 
we use have two, three or four degenerate minima, and all are non-negative functions of the 
field p. 
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Figure 1. Plot of the potential with four degenerate minima (3.1). 


3.1 Four degenerate minima 

The shape of the potential can in this case be parameterised as 

V(<p) = A^(<p^-ay(<p^-br, 


(3.1) 


where the constants are 0 < a < 6, A > 0. The potential (3.1) has four degenerate minima 
Ipi = —b, Tp2 = —a, ^3 = a, and ^4 = b. Following ref. [39], we use 

-1 + V3 ^ 1 + V3 

a= - - -, b =—-—. (3.2) 

We also set A = 1 in the numerical calculations, which amounts to measuring the potential 
V{ip) in units of A^, while x and t — in units of A“^. Fig. 1 shows the potential (3.1) 
corresponding to the above choice of parameters. For this potential, Eq. (2.5) gives: 

dip 


y/2\x = J 


(3.3) 


v /((^2_ o 2)2((^2_52)2- 

Looking at Fig. 1, one can deduce that the kinks that connect the neighbouring vacua 
belong to either of the three topological sectors (—6,—a), (—a, a), or (a,b). For example, 
a static BPS-saturated solution of Eq. (3.3) interpolating between the vacua (p = —a at 
X —)■ —00 and y? = a at X — 7 - +00 belongs to the sector (—a, a). At the same time, there 
is also the corresponding antikink that connects the same vacua, but (p = a at x ^ —00 
and (p = —a at x —)• +00. Formally, it belongs to the topological sector (a,—a); the 
distinction between kinks and antikinks is, however, just a matter of convention, and we 
will not distinguish between the sectors (—a, a) and (a, —a) and will drop the prefix “anti”, 
unless the opposite is needed in order to avoid confusion or to explicitly identify the field 
conhguration in question. 

The kinks corresponding to the three topological sectors above can be obtained from 
Eq. (3.3) as implicit functions of x [35, 39]. Below, we examine their excitation spectra. 

3.1.1 Topological sector (—a, a) 

This kink is constrained by \ip\ < a, which allows Eq. (3.3) to be rewritten as 

dip 


\/2Ax = J 


(a^ — (/?2)(52 _ 
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The corresponding implicit solution is [35, 39]: 


_ a + ip /6 - <y9\ 

a - if \b + ipJ 


(3.4) 


where p = 2\/2Aa(6^ — a^). This equation can be solved for ip{x) numerically, and the plot 
of this kink is shown in Fig. 2 (left panel). The corresponding antikink can be obtained 



U 



Figure 2. Left panel — the solution of Eq. (3.4) that connects the vacua —a and a. Right panel 
— the potential U{x) corresponding to the kink shown in the left panel. The dashed line shows the 
value of ujf. 


from Eq. (3.4) by substituting x —>■ —x: 

^ a + ip f b-ip Y'/^ 
a-ip\b + ip) 

Using Eq. (2.2) or Eq. (2.9) gives the energy (mass) of the static kink: 


(3.5) 


M{-a,a) = - O?). 

We performed a numerical search of excitations of the kink (3.4) lying in the discrete 
part of the excitation spectrum. Fig. 2 (right panel) shows the corresponding potential 
U (x) that enters the Schrodinger eigenvalue problem. This problem was solved using the 
standard methods, namely, integrating Eq. (2.13), with the known asymptotic behaviour of 
its solutions at x —)> ±oo, starting at a large negative x = xi (the “left” solution) and a large 
positive X = Xr (the “right” solution). The two solutions were then matched at some point 
X close to the spatial origin. The specific choice of x is not very important, for instance, 
one could take x = 0, however, it is convenient to take a small offset from the zero value 
— this helps to avoid technical issues when U{x) is an even function of x (in this case, 
excitation profiles ipix) are either even or odd functions of x, and the latter ones have a 
node at X = 0). We selected those values of uj at which the Wronskian of the “left” and the 
“right” solution, calculated at the matching point, turns to zero. 

We found two eigenvalues, Wg = —2 x 10“®, and iof = 2.70491. The former value is 
just the translational mode, whose exact energy is cjq = 0; the deviation of our numerical 
result from zero thus provides an estimate of accuracy. The latter value corresponds to 
the vibrational excitation, whose existence is non-trivial and reflects itself in resonance 
phenomena occurring in kink-antikink collisions in this topological sector, see Sec. 4. 
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3.1.2 Topological sector (—6,—a) 

In this sector a < |y5| < b, which turns (3.3) into 

Taking the integral results in 

if + a \b- if J 


(3.6) 


where p = 2\/2Aa(6^ — a^); note that the kink that connects the vacua a and b can be 
directly obtained from this expression. The static kink (3.6) has the mass 

^(-b-a) = A(6 — a)^(a^ + 3a6 + 6^); 

its profile is shown in Fig. 3 (left panel). 


V U 


- 0.4 
- 0.6 
- 0.8 
- 1.0 

_ -y 

- 1.4 

Figure 3. Left panel — the kink (3.6) that connects the vacua —b and —a. Right panel — the 
potential U(x) corresponding to the kink shown in the left panel. 

As opposed to the sector (—a, a), we found the kink (3.6) of the sector (—6, —o) (which 
applies to that of the sector (6, a), too) to only have the trivial translational excitation, 
the numerical value being Wq = — 7 x 10“®. The corresponding potential entering the 
Schrodinger equation is shown in the right panel of Fig. 3. It has to be noted that, even 
though the kinks considered in this subsection do not have a vibrational excitation, a static 
kink and the corresponding static antikink located close to each other can still have a 
vibrational mode. Notice that even though such a configuration is not a solution of the 
equation of motion, due to the nonlinear character of the field interaction, the introduced 
error typically falls off exponentially with increasing separation between the two solitons, 
which is why such Ansatze are routinely employed. The existence of such a “collective” 
vibrational mode can manifest itself as resonance phenomena in kink-antikink collisions. 
Realisations of this mechanism have been studied, for instance, in refs. [40, 41]. 

3.2 Three degenerate minima 

This situation corresponds to the potential taking the form 

V(<p) = AV(^^-ayy^^ + by, 
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(3.7) 



















Figure 4. Plot of the potential with three degenerate minima (3.7). 


where the constants satisfy a > 0, 6 > 0, A > 0. This potential has three degenerate 
minima: = —a, Tp 2 = 0, and = a. The corresponding topological sectors are (—a, 0) 

and (0, a). Following ref. [39], our choice of parameters is 


a = 


3 

4’ 


6 = 1 . 


Fig. 4 shows the plot of the resulting potential. Obviously, the sectors (—a, 0) and (0, a) 
are related by a change of the sign of (p{t,x), hence only one of the two sectors has to be 
studied. Eq. (2.5) with the potential (3.7) results in the following implicit expression for 
the kink that connects lp 2 — ^ ^3 = 


/ \/62 + ip'^ -b\^ ^ ^ 

\/P~+~a2 - ^/P~T^/ + bj 


where p = 2\/2Aa^\/P~+~P, with the mass of this kink being 

M(0,a) = (2(6^ + a2)5/2 _ + sP)) . 

Fig. 5 shows a plot of the kink (3.8) in the left panel and of the corresponding Schrodinger 
potential in the right panel. Our study of the excitation spectrum of this kink found only 
the translational mode, with the numerical result for the eigenvalue Wq = 4 x 10“^. 


3.3 Two degenerate minima 

The self-interaction potential of the model can in this case be written as 

V{p) = \\p^-a^f{p^ + h^)\ (3.9) 


where a > 0, 6 > 0, A > 0. The two degenerate minima are Tp-^ = —a and Tp 2 = a, and 
there is only one kink that connects the points —a and a, and its antikink counterpart. The 
parameters we use are [39]: 

4 
5’ 
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a = 


6 = 1 . 
















u 






Figure 5. Left panel — the kink (3.8) that connects the vacua 0 and a. Right panel — the 
potential U(x) corresponding to the kink in the left panel. 


V 



Figure 6. Plot of the potential with two degenerate minima (3.9). 


The potential (3.9) with these parameters is plotted in Fig. 6. 

The static kink belonging to the topological sector (—a, a) is, as usual, obtained from 
the equation of motion (2.5): 





\a-<pj 


(3.10) 


where ^ 


2y/2\a{b‘^ + a^). The mass of this kink is 


M{-a,a) = 

Fig. 7 shows the kink (3.10) in the left panel, and the corresponding quantum-mechanical 
potential that dehnes small excitations of this kink in the right panel. 

The kink of this specihc realisation of the model yields the richest set of eigenmodes: 
besides the translational mode (with the numerical value of Wq = 5 • 10“®), we found three 
vibrational excitations with energies = 4.27575, wl = 10.1893, and wl = 13.6095. 

A few more words are in order regarding the numerical accuracy. As stated above, 
the deviation of the numerical result for uJq from zero can serve as an estimate of the 
numerical error of the energy values. On the other hand, regarding the precision of the 
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Figure 7. Left panel — the kink (3.10) that connects the vacua —a and a. Right panel — the 
potential U{x) that corresponds to this kink. The dashed lines show the values of a;|, and a;|, 
ordered from bottom to top. 

excitation profiles 'ipix), one can, for instance, check that the profiles corresponding to 
different eigenvalues (where there are excitations other than the translational mode) are 
orthogonal. The results of this check are given in Table 1. 

Table 1. Eigenvalues and orthogonality checks (where applicable) of the kinks considered in Sec. 3; 
ipiix) are normalised to unity. 



Kink 

XI 

Xr 

Eigenvalue 

Orthogonality check 

(3.4) 

-9 

9 

~ -2 X 10“® 

— 

ujj ~ 2.70491 

(V’O, V’l) ^ -1 X lO”"^ 

(3.6) 

-9 

9 

CJg ~ — 7 X 10“® 

— 

(3.8) 

-42 

13 

Wq — 4 X 10“® 

— 

(3.10) 

-10 

10 

cjq ~ 5 X 10“® 

— 

ujj ~ 4.27575 

(^/^0,V’i) - 1 X 10"^° 

~ 10.1893 

(V'o,V’2) ^ 4 X 10"® 
(V^l,'02) ^ 2 X 10“® 

ul ~ 13.6095 

(V’o, V’s) ^ -4 X 10“® 
(V^i,V’3)-2x 10-5 

(V^2,V’3)-2x 10-5 


4 Kink-kink collisions 

The existence of a vibrational excitation of a kink means that, in certain conditions, kinetic 
energy of the moving kink can be transferred to and stored in the form of small oscillations of 
the kink’s profile. This is known to lead to interesting phenomena in kink-kink collisions [3, 
24, 40-44]. 

As an illustrative example, we performed a numerical calculation of the kink-antikink 
scattering in one of the topological sectors of the variant of the y?® theory with four degener¬ 
ate minima, i.e., with the field self-interaction given by Eq. (3.1). We studied a configuration 
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constructed from the kink (3.4) and the corresponding antikink (3.5). Namely, we set the 
initial field profile to 


(/9(f,x) 


f{—a,a) 


X + Xq- Vint 


1 — vt 


+ ¥^(a,-a) 


X - Xo + Vint 


1 — Vn 


- a, 


(4.1) 


where <^(_a,a)(®) the static kink and the static antikink. This setup 

corresponds to the kink and the antikink, separated at t = 0 by 2xo, and moving towards 
each other with the velocities ±Uin in the laboratory frame. As the separation between the 
kinks increases, the overlap between them becomes exponentially small and they stop being 
affected by each other — this happens in the asymptotic regime of the collision. In practice, 
the initial separation 2xo has to be much larger than the typical width of the kink; in our 
calculation, we use 2xo = 25. 

We solved the equation of motion using the standard explicit finite difference scheme, 

— ^2 ’ ~ 1^2 ’ 

where r and h are, respectively, the time and space grid spacings, and {k,j) number the 
corresponding coordinates of the grid points, {tk,Xj). The initial conditions follow from 
Eq. (4.1), and the presented results correspond to r = 0.002 and h = 0.01. The infinite 
space domain was truncated, so that —l<x<l, with I = 100, whereas the time varied 
in the range 0 < t < 900. As dictated by the dependence domains, the calculation started 
from a much larger interval of the x axis at t = 0 (namely, —4650 < x < 4650). To check 
our numerical results, we tested the conservation of energy as the time evolution progressed, 
taking into account energy flux through the endpoints of the space interval: 


1/a^y 1 
2 \ dt J 2 





dx 


dip dp 

dt dx 



2M 


where tc is the current moment of time (the integrand in the first integral is evaluated at 
t = tc). 

Our numerical simulations showed several different scattering regimes can realise, de¬ 
pending on the initial collision speed. We found out that there is a critical speed Vcr — 0.3160 
such that at Vin > Vcr the kinks bounce off each other and escape, as illustrated in Fig. 8. 
The kink-antikink scattering is only approximately elastic, i.e., the asymptotic escape ve¬ 
locity V{ < Vin- This feature can be seen in Fig. 8 as well; the typical values are, e.g., 
Vf = 0.23 at Vin = 0.40, Vi = 0.36 at Vin = 0.50, and Vf = 0.47 at Xm = 0.60. 

As opposed to the previous regime, collision velocities Vin < Vcr result in the formation 
of a long-lived bound state of the two kinks, the bion, see Fig. 9 (as observed long ago, for 
instance, in the model [45]). This, however, does not hold for all Xjn < Vcr', there are 
intervals of the collision velocities — the escape windows — where kinks still escape to the 
spatial infinity albeit only after having collided two, three, or more times. Similar processes 
have been observed in other field models, see, e.g., [3] for review. This phenomenon has 
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Figure 8. (Inelastic) kink-kink scattering at tijn = 0.5000. Left panel — space-time evolution of 
the system. Right panel — ip{t, 0), the time profile of the field at a; = 0. 




Figure 9. Bion formation at wjn = 0.2740. Left panel — space-time evolution. Right panel — 
the profile of (p{t,0). 


been explained by resonant energy exchange between the translational mode and a localized 
excitation mode of the kink. The value of the corresponding resonance frequency ujr in some 
models coincides with the kink vibrational mode frequency uji , whereas in other field models 
ujR could significantly deviate from oji [41, 42], or even belong to the continuous part of the 
excitation spectrum [43]. During the first collision, the kinetic energy of the kinks can partly 
be transferred to the mode IjOr. After that, the kinks bounce but cannot escape, and hence 
they stop at some point and then collide again. However, the second collision may result 
in transfer of the energy of the mode ojr back to the kinetic energy. This can happen if the 
time Ti 2 between the collisions and the frequency ojr are in a certain resonance relation 
(see below), and, as a result, the kinks can escape after the second collision. An example of 
such a two-bounce is shown in Fig. 10. Note that this regime is also inelastic — the kinks 
lose a part of their initial kinetic energy; this can be seen in the figures as well. We found 
around twenty two-bounce escape windows, as well as a number of three- and four-bounce 
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Figure 10. Two-bounce at Ujn = 0.2868. Left panel — space-time evolution. Right panel — 
the profile of 0). 


escape windows (where the kinks escape to inhnity after three and four collisions). Fig. 11 
shows the positions of these (/j® model escape windows and demonstrates that they exhibit 


N 

4 

3 

2 

1 


0.26 0.27 0.28 0.29 0.30 0.31 


Vin 



Figure 11. Left panel — Quasi-fractal structure of the escape windows; here N is the number 
of collisions before escape. Right panel — locations of two-bounce escape windows versus n, the 
number of small oscillations (see text for the definition). 


features similar to observed, e.g., in the model [3], namely, the quasi-fractal structure 
and the concentration of resonances near the critical velocity Vcr- We also analysed the 
positions of the two-bounce escape windows, using phenomenological relations that have 
been applied to resonance kink-kink collisions in other models, cf. ref. [3]. The main results 
of this analysis are shown in Fig. 12 and in Table 2. 

The hrst of the relations in question connect the resonance frequency ur with the 
number n of small oscillations of the held at x = 0 between the two collisions (so that, e.g., 
the two-bounce shown in Fig. 10 corresponds to n = 2): 

^rTi 2 = 27r(n -|- 2) -|- h = 2'nh + 6, (4-2) 

where (5 is a constant phase shift. Note that n can also be viewed as the number of the 
respective escape window in the sequence shown in Fig. 11 (implying that the windows with 
numbers 14, 15, 16, and 18 are missing — we have not observed them in our numerical 
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simulations). The shift of n by two in Eq. (4.2) is dictated by our choice that 5 has to be 
between 0 and 27r. Fig. 12 shows T 12 as a function of h. As one can see from this figure, 
the linear dependence of Eq. (4.2) fits nicely the numerical results, with the values of the 
resonance parameters resulting from the fit being ur = 1.618 and 5 = 3.403. This value of 
Ur is slightly less than the frequency of the kink’s vibrational excitation, ui = 1.644. This 
can be attributed to the interaction between the two kinks: it continuously distorts the 
excitation spectrum of the kink-kink system during the collision, which (inter alia) can lead 
to a deviation of ur from ui [41]. As mentioned above, the collective excitation spectrum 
of the kink-kink system can have a vibrational mode even when either of the solitary kinks 
does not. This feature gives rise to resonance phenomena such as escape windows and quasi¬ 
resonances in models such as where the kinks only have the translational excitation [40] . 


Tn 



Figure 12. Time between the first two collisions ri 2 as a function of h for two-bounce escape 
windows. Dots show results of our numerical calculations, the solid line is the corresponding linear 
fit of Eq. (4.2). 


Table 2 collects information we have obtained on two-bounce escape windows. It lists 
the initial collision velocities Uin (taken at the midpoint of the corresponding window), the 
final escape velocities Uf, and the values of T 12 . We also applied the phenomenological 
analysis of ref. [43], calculating 

/3 = Ti2y/vl-vf^ (4.3) 

for each of the escape windows, with the mean value being ^ ~ 2.569. The latter quantity 
was used to predict the locations of escape windows: 


theor 

^in 


(n) 


- '^CI 




R 


2ucr(27rh -|- 5)2 ’ 


(4.4) 


where is the predicted escape window initial velocity, assuming small deviation of 

the latter from the critical velocity, Uin — Vct- Table 2 shows that Eq. (4.4) does not 
provide accurate predictions for Uin — theoretical values at large values of n can differ from 
the corresponding calculated values by an amount larger than the distance between two 
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adjacent escape windows. This is at least partly due to poor accuracy of the determination 
of (3, which has a scatter of about 15%. We have thus not been able to use Eq. (4.4) 
to pinpoint the location of the escape windows with n = 14,15,16,18 that had not been 
observed directly in the numerical calculations. 


Table 2. Two-bounce escape windows. The values of Uin are given at the midpoint of the corre¬ 
sponding escape window. 


n 

Vin 

V{ 

Ti2 

(3 

..theor 

^in 

1 

0.2639 

0.158 

13.446 

2.337 

0.2608 

2 

0.2868 

0.207 

17.422 

2.311 

0.2824 

3 

0.2967 

0.229 

21.424 

2.330 

0.2934 

4 

0.3019 

0.239 

25.386 

2.370 

0.2998 

5 

0.3051 

0.237 

29.256 

2.407 

0.3038 

6 

0.3073 

0.247 

33.280 

2.451 

0.3065 

7 

0.3087 

0.214 

36.940 

2.495 

0.3084 

8 

0.3099 

0.253 

41.180 

2.545 

0.3098 

9 

0.3108 

0.079 

45.528 

2.599 

0.3108 

10 

0.3113 

0.222 

48.678 

2.643 

0.3116 

11 

0.3119 

0.253 

52.892 

2.684 

0.3122 

12 

0.3123 

0.249 

56.624 

2.730 

0.3127 

13 

0.3127 

0.124 

60.058 

2.736 

0.3131 

17 

0.3136 

0.041 

75.466 

2.934 

0.3142 

19 

0.3140 

0.221 

83.542 

2.965 

0.3145 


As mentioned above, in our calculations we also observed three-bounces and four- 
bounces — in these cases the kinetic energy is (partly) restored only after three and four 
collisions, respectively. Typical space-time pictures of these processes are shown in Figs. 13 
and 14. Another resonance feature, observed previously in kink-kink collisions in the double 
sine-Gordon [41, 42] as well as in the modified sine-Gordon model [43], is quasi-resonances. 
They occur when the energy that is pumped back into the translational mode during the 
second collision is not enough for the two kinks to escape. The two kinks merely come some 
large but finite distance apart and then collide again, eventually merging into the bion. This 
situation shows as a peak on the plot of T 23 , the time between the second and the third 
collisions, as a function of the initial velocity Uin. In some models, quasi-resonances have 
been shown to replace some of the escape windows [41], which is not unexpected as both 
these phenomena are due to the resonant energy exchange. Our calculation has not shown 
any clear signs of quasi-resonances occurring in the considered variant of the y?® model. 
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Figure 13. Three-bounce at V{n = 0.3012. Left panel — space-time evolution. Right panel — 
the profile of 0). 




Figure 14. Four-bounce at Ujn = 0.2568. Left panel — space-time evolution. Right panel — 
the profile of 

5 Summary and discussion 

We studied excitation spectra of the f/?® field model, choosing three different sets of the 
model’s self-interaction potential. We found that some of the kinks that arise in those 
variants of the model have vibrational excitation modes, which points at the possibility 
for resonance phenomena to occur in low-energy kink-kink scattering. To illustrate that 
relation between the excitation of a solitary kink and the kink-kink scattering, we performed 
numerical modelling of the latter process, using one of the considered kinks that has a 
vibrational excitation, and the corresponding antikink. We showed that there are two 
collision regimes, depending on the initial velocity of the two colliding kinks. Namely, the 
kinks bounce off each other inelastically at Uin > Ucr, whereas v\n < Vcr results in the 
formation of a bound state of two kinks, the bion. Furthermore, we demonstrated the 
existence of escape windows — intervals of initial velocities in the domain Uin < Vcr where, 
even though these velocities are below the critical speed, the kinks escape to infinities 
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after two, three, four etc. collisions. Our analysis of two-bounce escape windows shows 
that their positions and structure are well described by the phenomenological resonance 
condition, with the corresponding resonance frequency being close to the frequency of the 
kink’s vibrational mode. This confirms the resonant energy exchange between the kink’s 
translational mode and its vibrational mode (or, more precisely, a localised excitation of 
the kink-kink system that is close to the vibrational mode of a solitary kink) as the driving 
mechanism that leads to the occurrence of escape windows in kink-kink collisions in that 
sector of the (/J® model. 

Note that the value of the critical speed Vcr — 0.3160 that results from our analysis of 
collisions between the kinks (—a, a) and (a, —a) corresponds to the particular choice of the 
model parameters a and b. One can expect that the critical speed will be a function of these 
parameters; this distinguishes the y?® model from such models as and where there are 
no free parameters in the model (or, more strictly, they can be eliminated by appropriate 
rescaling of the field and the space-time coordinates). The critical speeds in these latter 
models are therefore fixed, namely, Vcr — 0.2598 in the (p‘^ model [3], whereas the model 
[44] has two different critical speeds (there are two different topological configurations of 
colliding kinks where the formation of a bion is possible and hence the critical speed can 
be defined), Vcv — 0.289 for the colliding kinks (0, —1) and (—1, 0) and Vcr — 0.045 for the 
kinks ( — 1,0) and (0,-1). In this connection, a more illustrative example is that of the 
modified sine-Gordon model [43], where the model potential has a single parameter r such 
that at r —)• 0 one recovers the regular sine-Gordon model (which is fully integrable and 
formally has Vcr = 0). The modified sine-Gordon model yields the following critical speeds: 
Vcr — 0.112 at r = 0.05, Vcr — 0.234 at r = 0.1, Vcr — 0.337 at r = —0.5. One can clearly 
see that the critical speed depends on the model parameter and that its value increases 
as the model moves away from the integrable limit. As far as the y?® model is concerned, 
it is probably safe to state that the critical velocity in the collisions that we studied can 
be made as low as that of the model. Indeed, one could consider a situation when b is 
chosen much larger than a. The relevant model potential, Eq. (3.1), would in this case be 
very close to that of the model as long as ~ a, leading to the dynamics of collisions 
between the kinks (—a, a) and (a, —a) being close to that of kink-antikink collisions in the 
ip'^ model. 

Finally, we would like to reflect on several issues that have been left out of the scope 
of this work but represent interest for future research. 

1. Some of the kinks of the model, corresponding to particular choices of the model 
parameters, can have power-law, rather than exponential, asymptotic at the spatial 
infinity [39]. Studying these kinks along the lines of our work would be a natural 
extension of this work. 

2. Resonant energy exchange in kink-kink collisions can also occur when the correspond¬ 
ing solitary kinks do not have vibrational excitations. For instance, two colliding kinks 
can have a collective vibrational mode, even though each of the two solitary waves 
does not. This has been shown to happen in the and the double-sine-Gordon 
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models [40, 41], This situation, apparently, can take place in some variants of the 
model. 

3. We think (in agreement with the authors of ref. [39]) that the collective coordinate 
method could be productively applied to studying kink-kink interaction in this model. 
It is worth noting that kinks with power-law asymptotic can apparently result in 
power-law asymptotic of the collective potential, which is a new feature that does not 
occur in other field models, such as, e.g., and </?®. 
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